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Abstract 

In the framework of the Inert Doublet Model, a very simple extension of the Standard Model, we study the 
production and propagation of antimatter in cosmic rays coming from annihilation of a scalar dark matter 
particle. We consider three benchmark candidates, all consistent with the WMAP cosmic abundance and 
existing direct detection experiments, and confront the predictions of the model with the recent PAMELA, 
ATIC and HESS data. For a light candidate, Mdm ^ 10 GeV, we argue that the positron and anti-proton 
fluxes may be large, but still consistent with expected backgrounds, unless there is an enhancement (boost 
factor) in the local density of dark matter. There is also a substantial anti-deuteron flux which might be 
observable by future experiments. For a candidate with Mjjm ^ 70 GeV, the contribution to e+ and p 
fluxes is much smaller than the expected backgrounds. Even if a boost factor is invoked to enhance the 
signals, the candidate is unable to explain the observed e+ and p excesses. Finally, for a heavy candidate, 
Mdm ^ lOTeV, it is possible to fit the PAMELA excess (but, unfortunately, not the ATIC one) provided 
there is a large enhancement, either in the local density of dark matter or through the Sommerfeld effect. 

1 Introduction 

The recent cosmological observations indicate that about 80% of matter in the universe is made of dark 
matter (DM) [1-3]. The most popular, particle physics, explanation is that DM is made of weakly interacting 
massive particles, or WIMPs [4,5], and the leading WIMP candidate is the neutralino, a supersymmetric 
spin 1/2 Majorana particle. However, spin one WIMPs, like in models with extra dimensions [6], and, to 
a lesser extent, scalar WIMPs, are considered as interesting challengers. An instance of the latter is the 
lightest stable scalar of the Inert Doublet Model (IDM), an extension of the Standard Model with two Higgs 
doublets and a discrete Z2 symmetry which is imposed as a simple way to prevent flavour changing neutral 
currents (FCNC) [7]. The scope and ambition of the IDM can not compete with those of the MSSM or 
models with extra dimensions, but the IDM is a simple model and, nevertheless, it does have an interesting 
phenomenology, as emphasised in [8,9], and e.g. [10]. Also, it encompasses some of the features of other 
models with scalar dark matter, like Minimal Dark Matter scalar candidates [11], or singlet scalars, like in 
hidden portal models [12,13]. 

The prospects for the direct and indirect detection of IDM dark matter, both through gamma rays from 
the Galactic Centre (GC), and through neutrinos from the Earth and Sun, have been addressed in [9,14-19]. 
In the present article, we study the production and propagation of positrons (e+) and anti-protons {p) 
that would result from its annihilation in the halo of the Galaxy The subject matter is timely, given 
the recent release of various observations of antimatter in cosmic rays, in particular from the Payload for 
Antimatter Exploration and Light-Nuclei Astrophysics (PAMELA) experiment, which has an excess in the 
positron fraction, e+/(e+ + e~) in the 10 — 80 GeV range [21] (but no excess in the p/p flux [22]), the 
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Advanced Thin Ionization Calorimeter (ATIC) experiment, which shows an excess of cosmic ray electrons 
(plus positrons) at energies of 300 — 800 GeV [23], and, finally, the High Energy Stereoscopic System (HESS) 
experiment, which has measured the flux of cosmic ray electrons (plus positrons) between 0.6 — 4 TeV [24]. 

These data, and their puzzling features, have generated a flurry of activity. As shown in various works 
[25-29], the PAMELA and ATIC anomalies may be of astrophysical origin, a natural source of positrons being 
nearby pulsars. Nevertheless, much of the activity has been focused on the possibility to interpret the excess 
of positrons as being due to annihilation or decay of dark matter particles in the galactic halo [30-40]. 

As we have already emphasised, the main advantage of the IDM is it simplicity. This is both a blessing and 
a curse. On one hand, it explains why the IDM should not pretend to fully explain the recent data. On the 
other hand, there is no analysis yet of antimatter production in the IDM, and we see this as an opportunity 
to fill a gap in the literature. 

The paper is organised as follow. The Inert Doublet Model is summarised in Section 2. Anti-proton 
and positron production and propagation are briefly presented in Section 3, together with a discussion on 
backgrounds expectations. Section 4 is devoted to the contribution to cosmic rays of annihilation of a IDM 
dark matter candidate. Cosmic rays features specific of the Inert Doublet Model and, in particular, their 
relevance for the results of PAMELA, ATIC and HESS experiments, are presented in Section 5. We give our 
conclusions in Section 6. 

2 The Inert Doublet Model 

The Inert Doublet Model is a two Higgs doublet model, noted Hi and H2, with an unbroken Z2 symmetry, 
such that 

Hi Hi and H2 -H2, 

and with all the other Standard Model particles even [8,9]. The potential of the IDM extension of the 
Standard Model may be written as 

V = i^l\H,\^ + fil\H2\^ + M\H,\* + \2\H2\^ + AalFiPlFap + X^IhIh^I' + ^ [{HlH2f + h.c] . (1) 

In this model. Hi contains the standard Brout-Englert-Higgs particle h (the Higgs for short in the sequel). 
The discrete symmetry, which prevents FCNC, is assumed to be unbroken, so that the H2 has no vacuum 
expectation value, and, consequently, the lightest Z2 odd particle is stable. The dark matter candidate may 
be either one of the neutral components of the doublet H2 = {H~^ , l/y/2{HQ + iAq))^ 

Ml = -2//? = 2Ali;^ 

Mho = M2 + (A3 + A4 + A5)^V2, 

Ml^ = /i^ + (A3 + A4 - A5)^V2 . (2) 

Depending on quartic couplings, either Hq or Aq can be the lightest particle. We choose Hq and, following [9], 
we define Xl = (A3 + A4 + A5)/2, which measure the trilinear coupling between the Higgs /i, and a pair of Hq. 
Choosing Aq instead would not change our conclusions. In our investigation of the model we choose /i2, A2, 
and the masses of scalar particles, including the mass of the Higgs, as input parameters. 

Experimental constraints on the IDM model from colliders have been discussed in [9], and further in [41], 
and [42]. The latter also discusses the prospect for discovery of the Aq and Hq at the LHC. In [41], the LEP I 
and II constraints on the neutralino are used to put constraints on the mass range of both Hq and Aq. These 
constraints are essentially summarised in the Figure 8 of that reference. The three benchmark dark matter 
candidates we will consider here, as well as the properties of the their companions Aq and H^, are all chosen 
to be consistent with LEP and WMAP data. 

Assuming that the Hq was in thermal equilibrium in the early universe, there are essentially three distinct 
mass ranges within which an IDM candidate has an abundance consistent with WMAP. In the sequel we 
refer to these ranges as the low mass (3 GeV ^ ^ 8 GeV) [17], the middle mass (40 GeV ^ ^ 80 
GeV) [9], and the high mass (500 GeV < Mi^, < 15 TeV) ranges [14]. 
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Candidates in the low mass range may annihilate only through the Higgs. The option of a light Hq 
is challenged, but not excluded, by current direct detection experiments (CDMS, Xenon and Cogent), and 
may even be compatible with the DAMA results [17]. Candidates in the middle mass range may annihilate 
through the Z (when co-annihilation with Aq is allowed) or through the Higgs. Some solutions in the middle 
mass range may also have a large annihilation cross section, through loop corrections, into a pair of gamma 
rays. This is relevant for indirect detection from annihilation at the GC, and such a gamma ray line would 
be easily observed by the GL AST/Fermi satellite [15]. Above ~ 80 GeV, annihilation into pairs 

is allowed with a large cross section, and the abundance from freeze-out is predicted to be well below the 
WMAP abundance [14]. At even higher masses, however, the annihilation cross section tends to decrease, 
av (X 1/M^^, and, consequently, the relic abundance increases, and, for some masses, may be consistent with 
WMAP observations. In this regime, an IDM dark matter candidate is analogous to that of Minimal Dark 
Matter [11]. However, in the IDM, there are more parameters to play with, and there is a whole range of DM 
candidates around 1 TeV with a cosmic abundance that is consistent with WMAP observations [14]. 



3 Antimatter in Cosmic Rays 
3.1 Anti-protons 

The anti-protons in cosmic rays are essentially by-products (secondaries). They are created by spallation, 
through collisions between protons, which are primary constituents of cosmic rays, and nuclei (essentially 
hydrogen) from the interstellar medium (ISM). The creation, destruction and propagation of anti-protons 
in the ISM and in the galactic magnetic field, may be described by a diffusion equation for the anti-proton 
density np. Solving this equation requires the distribution of matter in the Galaxy. Its shape is modelled by 
two cylindrical slabs, one representing the position of the ISM gas, and the other one the halo of dark matter 
(see e.g. [43]). We take these slabs to have radial extension Rg and heights 2hg (galactic disk) and 2hh (halo) 
respectively. Furthermore, the flux of cosmic rays is supposed to be time independent over the relevant time 
scale, so that only a stationary solution of the diffusion equation is required. There exist various models 
of anti-proton propagation, and they are all more or less in agreement with each others and with the data, 
including PAMELA [43-46]. 

In the present article, to compute the flux of anti-protons at the top of the Earth atmosphere, we have 
used the propagation model implemented in the DarkSUSY package [47]. We have interfaced DarkSUSY 
with microMEGAs2.2 [48,49], a versatile package that allows to compute the dark matter abundance (and 
the relevant branching ratios) of any WIMP dark matter candidate in models with a Z2 parity. 

In the DarkSUSY code, the stationary anti-proton densities are taken to satisfy the transport equation 

-^{E,x)=V ' [D{R,x) Vnpj - V '{u{x) np)-p{E,x) np^Q{E,x) = 0. (3) 

The description of propagation of anti-protons is accomplished by the first three terms. The first one describes 
diffusion in the galactic magnetic field. In each part of the Galaxy the diffusion coefficient D is assumed to 
be isotropic 

D{X) = D{Z) = DgO{hg - \Z\) + DhO{\z\ - hg), (4) 

and to depend only on the magnetic rigidity parameter R of the anti-protons (in Giga Volts or GV), 

A(it^) = A'(l + £j , (5) 

where I = g^h and Rq ^ 1 GV. The second term represents large scale convective motion, with a velocity field 
u{x). This term is introduced to model the impact of the wind of cosmic rays, that blows away from the disk, 
on the motion of anti-protons. The third term, finally, models the loss of anti-protons due to collisions with 
the ISM, mostly hydrogen, with 

p(E,x)=nH{x) v{E) a{E), (6) 

where a is the inelastic cross section, and uh is the number density of hydrogen in the galaxy, which is 
assumed to be of the form 

nH{x) = nniz) = nHgO{hg - \z\) + nHhO{\z\ - hg) . (7) 
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As for boundary conditions, it is assumed that the density of cosmic rays, and thus of anti-protons, is negligible 
at the boundary of the Galaxy, i.e. 

np{Rg, z) = np{r, hh) = np{r, -hh) = . 

Finally, the last term on the r.h.s. in ([3| is the source of anti-protons with energy E. Possible sources 
include secondary anti-protons produced by spallation and, possibly, dark matter annihilation or decay in 
the galactic halo. To discriminate between standard, astrophysical sources of anti-protons, and more exotic 
possibilities, like dark matter, it is clearly necessary to have a good handle on the expected background 
signal. For low energies, the transport models sketched above give a flux of anti-protons that is in good 
agreement with observations, in particular if p — He collisions and energy loss effects during propagation are 
taken into account [46]. However, the uncertainties on the various parameters (transport coefficients, nuclear 
processes,...) imply that the background flux at high energies may not be uniquely determined from the low 
energy data [50]. Nevertheless, for the sake of the argument, and like in most analysis of antimatter in cosmic 
rays, we assumed that the background is known. More specifically, we adopt the flux used in the recent 
analysis of [50] 

wri^TWmimW^ [Gev-..,.-. (8) 

where Tp = Ep — rrip and t = Tp/1 GeV. 

The anti-proton data are generally presented in terms of the ratio of the anti-proton to proton fluxes, 

p ^ d(l)/dTp 
p d(j)/dTp 

Thus we need the flux of protons. For the background, we use the spectrum [51] 

g: = 8+LlA"+o.8,..»» [GeV-.-m-^.-.s,-.]. (10) 

The ratio p/p observed is small, typically (9(10~^), thus it is sufficient to only consider the background proton 
flux in the ratio (|9|. 



3.2 Positrons 

Like anti-protons, positrons in cosmic rays are supposed to be secondaries. Positrons may be created in the 
decay of pions and kaons, which themselves are generated through spallation. Like anti-protons, the motion of 
positrons in the galactic magnetic field is diffusive. They may also loose energy through synchrotron radiation 
in the interstellar magnetic field, and through inverse Compton scattering on diffuse starlight or on the cosmic 
microwave background. 

We follow here the classic treatment of [52] . The positron number density per unit energy, dn^^ / dE^^ , is 
given by the stationary solutions of a transport equation 



d dnQ+ 
didE^ 



= V- 



d 



(11) 



where K{Eq+^x) is the diffusion constant, b{E^+^x) is the rate of energy loss, and Q(£^g+,x) is the source 
term (in cm~^ The relation between the flux and the density number is 



d^ 



(3c dUe 



dE^+ 47r dE^+ 



(12) 



It is generally assumed that the diffusion constant is independent of position within the diffusion zone, and 
its energy dependence is taken to be given by 



K{E) = Ko 



C 



IGeV/ 



(13) 
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where Kq = 3 • lO^^cm^s ^, C = 3" implements a plateau in the diffusion constant below 3 GeV, and the 
spectral index is = 0.6. [52]. The energy loss rate depends on the positron energy through 




where te = lO^^s is the characteristic time for energy loss. The diffusion zone is modelled in this case by a 
slab of thickness 2L with L = 3 kpc. The diffusion equation is solved by considering a leaky box model, with 
free escape boundary conditions, i.e. the cosmic ray density is set to zero on the boundary of the slab. 



Again, to confront the contribution of dark matter annihilation to positrons in cosmic rays, we need a good 
handle on the expected e+ background. An extensive discussion of the many uncertainties on e+ production 
by spallation and their subsequent propagation in the interstellar medium may be found in [53]. A further 
issue is that, to interpret the data, it is also necessary to know the electron flux. This is because the positron 
data are generally presented in term of the positron fraction 

e+ 



dE 


e+ 


d^ 
dE 


I d^ 

e+ dE 


e~ 



This ratio permits to factor out the, energy dependent, acceptance of the detectors (supposed to be the same 
for electrons and positrons), and, to some extent, the effect of solar modulation, to be discussed below. 

Unfortunately, the electron flux spectrum is poorly known. First, it is difficult to simulate the e~ flux 
in cosmic rays from first principles, simply because there are many astrophysical sources of electrons in the 
Galaxy, starting with our Sun. Second, as of today, the experimental situation is not very clear, as one may 
appreciate from the compilation of existing data presented in [54] , which show a wide spread in spectral index 
and even flux. Things will improve in the near future, with the advent of new data, including those on the 
spectra to be released by PAMELA. In the meantime, we have to confront the fact that the electron flux 
uncertainties have a non- negligible impact on interpreting the positron flux data, as discussed in [53]. In the 
present article, were we focus on particle physics predictions, we adopt a common practice and take a fixed 
background electron and positron fluxes. Concretely, we follow [52], and use the following background flux of 
positrons and electrons. 



dE 

dE 

M 
dE 



16 f~^-^ 

[GeV-i- cm-2. g-^- sr-^] (15) 



e , prim. 



1 + 11 eO-9 + 3.2 e2 i5 
0.70 eO-^ 



1 + 110 ei-5+600 e2-9 + 580 e^ '^ 
4.5 



[GeV-i- cm-2. s-^- sr-^] (16) 



e+, 



1 + 650 e2-3 + 1500 e4-2 



[GeV-i- cm-2. s-^- sr-^] (17) 



where e = (£^/lGeV). As for we compute the flux of positrons from IDM dark matter annihilation using 
DarkSUSY. 



3.3 Solar modulation 

Before closing this section, we mention one further complication. Low energy {E ^ 1 — 10 GeV) cosmic 
rays (both e+ and p) are strongly affected when they enter the Solar System, because they may loose energy 
by interacting with the solar wind and the solar magnetic field. This is generically called solar modulation. 
The standard approach assumes that solar modulation is independent of the sign of the charge, and thus 
the same for electrons and positrons (for discussion see [52,55]). Using the, so-called, Gleeson and Axford 
analytical force-field approximation [56,57], the flux at the Earth d^^/dE^ can be deduced from the flux at 
the heliosphere boundary, or interstellar (IS) flux d^js /dEjs, by 

dE^ pjsdEjs' ^ ^ 
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where the energy at the hehosphere boundary is given by 



Eis=E^ + \Ze\cl>. 



(19) 



Here and pis are the momenta at the Earth and at the hehosphere boundary respectively, e is the absolute 
value of the electron charge, and Z is the charge of the particle, in unit of e. In this approximation, the solar 
modulation is completely determined by the rigidity parameter ^, i.e. the particle momentum divided by its 
charge, expressed in volts, [(j)] = V. Since this solar modulation effect is the same for positive and negative 
charge particles, it does not affect the e+ fraction and the p/p ratio. 

However, the recent data depart from this simple picture, and reveal some effects which may be attributed 
to a sign dependent solar modulation effect. In particular, below 10 GeV, the PAMELA experiment has 
measured a substantially smaller positron fraction than the previous experiments (essentially CAPRICE and 
HEAT [58-60]). Similarly, over the years, BESS has measured variations in the p/p ratio. This sign dependent 
solar modulation effect is believed to be related to inversion of the solar magnetic field polarity, that take 
place with a periodicity of 22 years, in phase with maximum of solar activity, that occurs with 11 years 
periodicity (for details, see for instance [61] and [62]). Still, for the time being, there seems to be no complete 
understanding of this extra solar modulation effect. This of course somehow limits the use we could make of 
low energy data (both on e+ and p), in particular to constraint models of dark matter. However we will try 
to argue otherwise. 

4 Dark Matter annihilation 

We now consider specifically annihilation, in the galactic halo, of an IDM dark matter candidate into Standard 
Model particles. The kind of annihilation products one may expect will depend on the mass range being 
considered. However, in the simplest version of the IDM, there is no direct annihilation into lepton-antilepton 
pair^ and a fortiori into positrons. Instead, positrons (and anti-protons) will results from decay of primary 
annihilation products, like VK+l^" or hh pairs. This basic feature implies that we may not expect an IDM 
candidate to give a very good fit to both PAMELA and ATIC observations, which instead favour dark matter 
candidates with a large branching ratio (BR) into leptons (see e.g. [30,37]). 

The annihilation of an IDM dark matter particle into positrons or anti-protons enters the diffusion equa- 
tions through the source term. 



where (crann'^)^o^7 ^iid are respectively the total (average) annihilation cross section times the relative 
velocity the BR into a state / and the fragmentation function of the final state / into positrons/anti-protons. 

We have implemented the IDM in the microMEGAs2.2 code, which gives the relic abundance of dark 
matter as a function of the model parameters. The code also gives the annihilation BR. Depending on the 
mass (and parameters), annihilation takes place dominantly either through the Higgs (low and middle mass 
ranges) or into W+I^" and Z pairs (middle and high mass ranges). The BR are shown in Figjl] for a Higgs 
mass Mh = 120 GeV, and two values of the bare mass parameter, /i2 = 40 and /i2 = 200 GeV. From Figjl] 
we see that above 80 GeV the Hq annihilates dominantly into weak gauge bosons. Once Mhq > M^, the 
annihilation into h pairs opens, with a rate which depends on the Hq — h coupling A^. 

The differences between the two panels of Fig. [l]are perhaps worth being explained. Below the threshold 
for weak gauge bosons production, the BR are the same in both figures. This is because, in this regime, the 
Hq may only annihilate through the Higgs and the BR depend only on the Yukawa couplings of the fermions. 
Differences appear above the threshold. In particular, in the right panel, there is a dip in the fermion BR and 
two dips in the hh BR. Both are related to the change of the coupling A^, which is positive for Mh^ > /i2- 
In the right pane|^ the first dip corresponds to the vanishing of Al at Mhq = /i2 = 200 GeV. The second dip 
in hh production is due to a destructive interference, which may occur for A/, > (see [14]). 

^An exception is a light WIMP candidate, Mhq ^ GeV, which may have a large branching ratio into fr pairs, but this range is 
not relevant for explaining the PAMELA and ATIC features. 

^In the left panel, Xl also vanishes, for Mhq = 40 GeV, but, below weak gauge bosons production, the BR do not depend on Al- 




(20) 



^ 10' 10^ 10' ^ 10' 10^ 10' 

M^o [GeV] M^o [GeV] 

Figure 1: Branching ratios of the Hq annihilation as a function of Mh^. The Higgs mass Mh is fixed at 120 
GeV. Left panel: fi2 = 40 GeV; Right panel: /i2 = 200 GeV. 



The rate of dark matter annihilation is proportional to the square of the number density n(r) = p{f) /Mdm^ 
which is not well known. Observations of rotation curves suggest a rather cored profile [63,64], with a flat 
behaviour at the centre, whereas numerical simulations predict more cuspy profiles in the innermost region of 
the galactic centre (see Kravtov et al. [64], Navarro- Frenk- White (NEW) [65] and Moore et al. [66], as well as 
the recent Via Lactea and Aquarius simulations [67,68]). In the present study, we focus on the popular NFW 
distribution, which is a sort of benchmark in the community. The NFW profile is parameterised as follows 



p{r) = po 



1 + jro/ao)'^ 
1 + {r/aor 



(21) 



where ro is the galacto-centric distance, po = 0.3 GeV/cm^ is the dark matter density in the solar neighbour- 
hood, and (ro [kpc], ao [kpc], a, /5, 7)= (8.0, 20, 1, 3, 1). 

On top of this smooth, averaged distribution, numerical simulations indicate that substructures or clumps 
of dark matter may survive in virialized systems [67-72], and it has been suggested that such clumps may 
enhanced dark matter annihilation rates. This effect is usually parameterised through a boost factor (BF) 

op. _ Iv P dumpy ^ ly ^ P dumpy /^^x 

ivd'-pi^oot,~ Ply ' ^ ^ 

where the integration is over the volume which contributes to the annihilation flux, typically a few kpc^. 
Recent developments tend to disfavour the possibility of large enhancement of the fluxes [73-75], and, fur- 
thermore, it has been shown that realistic boost factors have an energy dependence [73,74]. Finally, it is 
expected to be different for positrons and anti-protons. Despite these potential caveats, we follow the standard 
practice and consider boost factors BF = (9(10 — 100) in this work. 

Another source of possible enhancement of the annihilation rates has recently received much attention. It 
has been shown that, in the presence of attractive long range interactions, the annihilation of non-relativistic 
dark matter particles may be enhanced in the limit of small relative velocities [76] . This so-called Sommerfeld 
aka Sakharov effect [77] may lead to spectacular enhancements, 0(10 — 10^), in particular near resonances. 
A general analysis of the Sommerfeld effect is quite complex [76] but, a rule of thumb is that it is more 
relevant for heavy WIMPS, M^m ^ TeV. More precisely, the large mass IDM case is very much analogous 
to the Minimal Dark Matter scenario, provided M^m ^ /o^w^ with W boson exchange giving effectively 
long range interactions with enhancements of (9(10^) [78]. For the low and middle mass range, however, the 
Sommerfeld effect is negligible. 
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Figure 2: Left panel: Positron fraction measured by balloon experiments and PAMELA (in red). Centre 
panel: Total electron and positron flux measured by balloon experiments, in particular ATIC, together with the 
HESS data. Right panel: p/p fraction. In all cases the red (blue) lines represent the expected background from 
Eq. (13-15) and Eq. (8-9) for leptons and baryons respectively without (with) a solar modulation of = 500 MV. 



5 The IDM vs the PAMELA, ATIC and HESS observations 

The PAMELA satellite has published its first scientific results on e+ and p in the Fah of 2008 [21]. The 
primary goal of PAMELA is the study of the antimatter component of moderate energy cosmic rays, with 
better resolution and, being a satellite and not a balloon experiment, with far better statistics than previous 
experiments. 

Concerning positrons, AMS 98 [79] and balloon experiments like CAPRICE 94 [80], CAPRICE 98 [58] 
and HEAT 00 [81], have taken data between 300 MeV and 40 GeV. PAMELA has been designed to extend 
this range from 50 MeV to 270 GeV. The present data published so far cover the range between about 2 and 
80 GeV (see Figj2j left), corresponding to about 10^ positrons identified [21]. 

The PAMELA collaboration has not yet released their observations of absolute flux of e~ and e+, but the 
PPB-BETS, ATIC and HESS collaborations have measured the flux of electrons in the energy range 30 GeV 
- 4TeV [24,82,83]. Unlike PAMELA, these experiments can not distinguish electrons from positrons, hence 
in principle their data (Figj2j middle) may also include positrons. The same figure also displays the recent 
HESS observations. 

The anti-proton flux has been analysed by IMAX 92 [84], BESS 95+97 [85], BESS 99 [86], BESS 02 [87], 
BESS Polar [88], CAPRICE 98 [89] and HEAT 00 [90] experiments, from 120 MeV to 40 GeV. PAMELA has 
been designed to improve this measurement for energies from 80 MeV to 190 GeV. The present published data 
give the p/p ratio up to 100 GeV, corresponding to about 10^ p identified (Fig|2j right [22,62]). Ultimately, 
PAMELA measurements will have a high statistics, with 10^ (10^) anti-proton (resp. positron) events after 
three years of operation. The data displayed in Figj2]are compared with the expected backgrounds, (8pQ) 



and ( T5p7 ). For anti-protons, there is a good agreement between expectations and observations, but the 



positron and electron data show significant deviations. 

First, compared with previous experiments, the PAMELA data show a depletion of the e+ fraction below 
5 — 6 GeV. This discrepancy is not fully understood yet but, as discussed above, it is assumed to be due to a 
sign dependent solar modulation effect, related to the periodic inversion of the solar magnetic field. Naively, 
we would expect an anti-correlation in time between the positron and anti-proton fluxes {i.e. minima of 
positron fraction corresponding to maxima of p/p ratio) but, this is not what data on anti-proton show. 
Fig. [2] However, simulations of propagation of low energy cosmic rays indicate that light and heavy particles 
are affected differently by sign dependent solar modulation effects [61]. This is manifest in the Fig. 2 of [61], 
which shows that both the e+ fraction and p/p ratio should be low at the time of PAMELA data taking, 
which is actually quite in agreement with observations (see also the talk by S. Ricciarini [62]). 

Second, the PAMELA data also show a steadily increase of the e+ fraction above ~ 10 GeV. The standard 
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Figure 3: Left panel : The thin curve represents the positron fraction generated by a fiducial flux of positron 
^{E) = 6.5 • 10^ E~'^-^ GeV"-*^ cm~^ s"-*^ sr"-*^. Right panel : Expected total electron plus positron flux in 
comparison with the ATIC measurements. In both panels, the thick plain lines represent the backgrounds, 
whereas (thick) dashed lines represents the signals (+ bgd). 

expectation, based on simulations of e+ production through spallation and propagation in the ISM, is that 
the e+ fraction should decrease instead of increasing with energy. This is manifest in the expected flux of 
Eq.([l5][l7|) which give 



(X E 



-0.15 



bckg 



The data (focusing on the three high energy points, see Fig. (|3| gives 



PAMELA 

This positron excess is the cause of the recent excitement around the PAMELA data, a possible explanation 
being dark matter particles annihilating in the vicinity of the solar system. Furthermore it is tempting to 
correlate the PAMELA excess to the electrons+positrons excess observed by ATIC at energies in the 300 — 800 
GeV range. This excess is also consistent with the HESS data. 

Actually, if we fit the PAMELA data by adding to the background an extra positron flux, with a simple 
power law spectrum, and an equivalent electron flux, the excess in PAMELA and ATIC data may be put in 
correspondence. This is suggested in Figjsj were we used ^e-,e+ E~'^-^. Our flux somewhat overshoots 
the ATIC data, but this could be ameliorated by lowering the expected background (plain blue curve in 
the right panel of Figj3|, something we were however reluctant to do. Obviously, the flux should also be 
cut-off at an energy around the ATIC peak. Altogether, such a spectrum may result from a source with an 
energy injection around the ATIC excess. These basic features are somehow consistent with the astrophysics 
explanation, which posits that the electron and positron excesses are due to cosmic rays produced by nearby 
pulsars [25-29]. Explaining both PAMELA and ATIC data from dark matter is in general more challenging 
but, as discussed extensively in the recent literature, a very simple, albeit ad hoc, culprit would be a heavy 
dark matter particle that annihilates or decays (in the case of a long-lived dark matter candidate) dominantly 
in lepton pairs. 

In the sequel of this work we investigate annihilation of an IDM dark matter into antimatter, and confront 
the model to the data. 



5.1 Low mass range 

In [17] it has been shown that a light IDM dark matter candidate {i.e. a light WIMP) may be consistent with 
both the WMAP and the DAM A/LIBRA results. DAM A/LIBRA is a direct detection experiment that has 
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Figure 4: A Low Mass candidate {Mdm = lOGeV). Thick Plain: Background, (Thick) Dashed: Signal (+ 
BckGrd) and (Thick) Dotted: Signal with boost factor = 10 (+ BckGrd). A solar modulation oi (j) — 500 MV 
has been applied. 




Figure 5: Left panel {Right panel) : Flux of positrons (anti-protons) for M^m — lOGeV. Thick Plain: Back- 
ground, (Thick) Dashed: Signal (+ BckGrd) and (Thick) Dotted: Signal with boost factor = 10 (+ BckGrd). 
The red (blue) curve corresponds to signals without (with) solar modulation (0 = 500 MV). 



reported evidence for an annual modulation of the nuclear recoils in their detector [91]. Taking into account 
the null results of all the other DM detection experiments (and the channelling effect on the threshold energy 
in DAM A), the DAM A/LIBRA signal may be due to elastic scattering of a dark matter from the galactic 
halo, with a dark matter particle mass in the range (see [92]) 

3 GeV< Mdm< 8 GeV. 



For these masses, both annihilation of the Hq and its scattering with a nucleus take place through the 
Higgs particle, with cross sections scaling like a ex Xj^/mf^. Hence, for a fixed DM mass, imposing WMAP 
abundance fixes all the relevant parameters of the model (modulo some residual uncertainty in the Higgs- 
nucleus couplings). That agreement with DAM A/LIBRA may be reached without further tuning is thus 

^The best fit to DAMA/LIBRA data gives a somewhat larger mass Mdm '^10 — 12 GeV, but these values are excluded by the 
CDMS and Xenon data [93,94]. 
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remarkable. This prediction is however not specific to the IDM. As emphasised in [17], any model with a 
scalar DM candidate coupled through the Higgs portal would give a similar fit to the data. 

In our discussion of positron and anti-proton signatures we consider as a benchmark, a light WIMP 
scalar candidate with mass M^m = 10 GeV and a si ~ 3 • 10~^^cm^. This candidate is somewhat within 
the boundary between the ranges allowed by DAM A/LIBRA and that excluded by all other experiments. 
However, the version of DarkSUSY code that we have adapted does not give the flux of positrons and anti- 
protons for a dark matter candidate lighter than 10 GeV. Computing the flux for the whole DAM A/LIBRA 
range would require to update the Pythia tables of DarkSUSY toward lower masses, but we have not done 
so. Still, we do not expect much difference between the flux produced by 10 GeV and 8 GeV dark matter 
candidates. Substantially lighter candidates, say with Mdm ^ 4 GeV, could give larger flux (albeit at lower 
energies), but these candidates are less favoured by the DAM A/LIBRA data. So, at this exploratory level, 
we limit ourself to a candidate with Mjjm = 10 GeV. 

A cosmic abundance consistent with WMAP is reached for ^ —0.2 for Mh = 120 GeV and Mjjm = 10 
GeV. Consequently, the flux of positrons and anti-protons is fixed, modulo the usual uncertainties regarding 
the distribution of dark matter in the halo or in the cosmic ray propagation parameters. The results are 
given in Figji] In all three figures, the dashed line corresponds to a flux computed for a standard NFW 
distribution, while the dotted line correspond to a signal boosted by a factor of 10. Given the low mass of 
the i^o, the signal is within the range of energies were both the positron and anti-proton fluxes are subject 
to solar modulation. This complicates the comparison between theoretical predictions and data. 

For both positrons and anti-protons, the contamination from dark matter annihilation is consequent only 
if we consider boosting the signal. Here we take BF = 10. Also the effect of the boost is more dramatic for 
anti-protons than for positrons. This is essentially because the expected background signal is much smaller for 
p/p than in the case of the positron fraction. Focusing on positrons, the excess from dark matter annihilation 
may be somewhat attenuated (hidden) by playing with the rigidity parameter (we have used (j) = 500 MeV 
in the plots), but doing so would also affect all low energy cosmic rays fluxes (including heavier nuclei) so 
we have refrained doing so. This still leaves open the possibility of hiding the dark matter signal with a sign 
dependent solar modulation effect (provided the boost factor is not too large) . It remains to be seen whether 
this can be done in a consistent way for all cosmic ray constituents, which differ in mass or charge. Moreover, 
dark matter annihilates into fermion-antifermion pairs, and so an excess would also be manifest in absolute 
fluxes. These are shown in Fig js] for positrons and anti-protons, together with existing data AMS 98 [79], 
CAPRICE 94 [80], CAPRICE 98 [95], BESS 95+97 [36], BESS 98 [96] and BESS Polar [88]. In the absence 
of a systematic analysis of (sign dependent) solar modulation in the light of recent data, for the time being 
it is probably fair to conclude that little may be said from positrons data alone. However, the current data 
on anti-protons allow to exclude BF larger than one. In the future, we expect constraints from positrons to 
remain much milder than those from anti-protons, because the DM contribution is typically smaller than the 
overall effect of solar modulation. We should emphasised that BF = 10 is already a large boost factor, and 
a fortiori so for anti-protons. Indeed, since anti-protons may travel much greater distances than positrons 
without losing energy, the anti-proton flux is expected to be less subject than the positrons to clumps in the 
dark matter distribution [75] . This means that the astrophysical BF is expected to be smaller for anti-protons 
than for positrons [74]. 

To conclude the discussion on a the light IDM dark matter candidate, we briefly comment on the production 
of anti-deuteron. The production of anti-deuteron by spallation is typically small and is predicted to fall off 
for kinetic energies below 1 GeV per nucleon. Given its large abundance, a light WIMP may give a substantial 
contribution to the flux of anti-deuteron at low energies (see for instance [97]). For the IDM candidate with 
Mjjm = 10 GeV, we obtain, using the DarkSUSY routines, an anti-deuteron flux at = 0.25GeV/n of 
9 • 10~^ (GeV/n s sr m^)~^ (for BF = 1), which is below the upper limit of 1.9 • 10~^ (GeV/n s sr m^)~^ 
set by the BESS experiment [98], but above the expected acceptance of the future AMS-02 and GAPS 
experiments [99,100], which are 4.5 • 10~^ (GeV/n s sr m^)~^ and 1.5 • 10~^ (GeV/n s sr m^)~^ respectively. 
Anti-deuteron data might turn out to give the strongest constraint on a light IDM dark matter candidate. 

5.2 Middle mass range 

In this subsection, we consider an IDM dark matter candidate in the mass range 

40 GeV< Mdm< 80 GeV. 
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E [GeV] E [GeV] T [GeV] 

Figure 6: Same as Fig. [i] but for a Middle Mass candidate (Mdm = 70 GeV, BF = 10^). 



In this range, both annihilation through the Higgs and co-annihilation through the Z are important to 
determine the relic abundance [9,14]. An illustration is given in Fig. [6j where the fluxes and fractions from 
the annihilation of a 70 GeV candidate, which has a cosmic abundance that is in agreement with WMAP. 

Since annihilations produce dominantly bb pairs, the positron spectrum emerging is very soft. Even if one 
considers a large boost factor to enhance the signal (BF = 10^ in the figure), the spectrum does not exhibit 
the steep increase in the positron fraction observed by PAMELA. There is a slight excess of positrons below 
10 GeV, but it is consistent with observations. Moreover this excess is within the energy range where solar 
modulation is effective, and may thus be easily hidden. 

Like in the case of a light WIMP, the excess of p is more significant than in positrons. Although the signal 
may also be mimicked by solar modulation for moderate boost factors, large boost factors (say BF ^ 10^) 
may be clearly excluded. Finally, there is a small excess around 10 GeV that perhaps could be probed with 
better statistics and better understanding of solar modulation effects. 

All these fiuxes have been computed for = 120 GeV. However it is of interest to envision larger Higgs 
masses, since this is one of the historical motivations for the IDM. In [9], it has been shown that Higgs as 
heavy as ~ 500 GeV may be consistent with LEP precision measurements. Although one expects that 
the annihilation cross section is small for such a large Higgs mass, it has been shown in [15] that the WMAP 
abundance may be reached using co-annihilation. In these scenarios, loop effects may be large and the dark 
matter candidate could have a large branching ratio into 7 pairs and into 7Z. Such gamma rays could be 
a promising signal for the GL AST/FERMI satellite. Although these properties require some amount of fine 
tuning, the prediction of significant gamma ray lines is specific to the IDM scalar dark matter candidate. Given 
that a Z may decay into lepton pairs, it is interesting to investigate whether loop effects may also enhance the 
production of positrons. To address this question, we have considered, as an instance, the benchmark Model 
I of [15], and have simply taken into account the positrons coming from Z decay. At one- loop, there is also a 
direct contribution with llj annihilation [36], but we have neglected it for our estimate. This approximation 
is a posteriori legitimated by the fact that we conclude from FigjT] that the signal in positrons from loop 
effects is too small to be observable, unless the boost factor is indeed very large, BF ~ 10^. Furthermore, 
although the shape of the spectrum is slightly harder for E ^ Mdm than it is at three level, loop effects may 
not reproduce the observed excesses and, at best, we could put a limit on the BF. 

5.3 High mass range 

For a heavier IDM dark matter candidate, annihilations into weak gauge bosons become dominant. Specifi- 
cally, in the IDM, we may consider dark matter candidates in the mass range 

500 GeV < Mdm ^15 TeV, 
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Figure 7: Same as Fig. [i] but for a Middle Mass candidate {Mdm = 70 GeV, BF = 10^) with one-loop 
contribution. 
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Figure 8: Same as Fig. [i] but for a High Mass candidate (Mdm = 1 TeV, BF = 5.10^ 
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Figure 9: Same as Fig. [i] but for a High Mass candidate (Mdm = 2 TeV, BF = 2.10^ 
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where the upper bound comes from unitarity hmit. In [14], it has been shown that the WMAP cosmic 
abundance requires the bare mass /i2 to be close to Mhq^ and thus, for most the parameter space, annihilations 
take place essentially into ly+I^" and Z pairs, just like for Minimal Dark Matter candidates [11]. The Figures 
([8|, ([9| and (10) illustrate the typical signals one may expect from IDM candidates with masses Mom = 
1, 2 and 10 TeV. Note that, in these cases, we have adopted an isothermal profile for the distribution of 
dark matter in the Galaxy, so as not to contradict constraints coming from synchrotron radiation emission 
by dark matter annihilation products [34]. The ATIC balloon experiment shows an excess between 300 and 
800 GeV which is consistent with the first data points of HESS. It is tempting to associate those signatures 
as being due to annihilation of a dark matter particle with a mass around TeV. For this to work requires a 
large boost factor BF ~ 10^, however, as for Minimal Dark matter, this boost may come from a combination 
of astrophysics and particle physics effects, the latter in the form of Sommerfeld enhancement. This being 
assumed, the strongest constraint comes again from p data from PAMELA, which require the candidate to 
have a small branching ratio into baryons (see e.g. [30]), or, focusing on the IDM, to be heavier than 10 TeV 
(FigjTol). Unfortunately, as it is well known, the fit to ATIC/HESS is then quite poor. 

One may envision improving the fit to the data by extending the IDM, for instance by adding more 
particles. This opens the possibility to have direct annihilation into charged lepton pair^ However, for 
a scalar DM candidate like in the IDM, these processes would in general be p-wave or chirality (ex mf) 
suppressed, and thus not very useful to explain the excesses. More sophisticated extensions may however be 
considered, as in [101], which do not have such a limitation. 



6 Conclusion 

We have confronted the IDM scalar dark matter candidate to the recent data on antimatter in cosmic rays. 
Not surprisingly, the fits to data are generically poor compared to those of the recent models that have been 
designed to explain the PAMELA and ATIC excesses. If one insists on explaining PAMELA and ATIC, the 
most promising possibility offered by the IDM is that of a very heavy DM candidate, with Mjjm ^ 10 TeV 
so as to evade the constraints on the anti-proton flux. In this regime, the predictions of the IDM are similar 
to those of a Minimal Dark Matter candidate, and the price to pay is the same, as both require a large boost 
factor to enhance the fluxes. The Sommerfeld effect, which is generically speaking relevant for a heavy dark 
matter candidate, may help, but still, some astrophysical boost factor would be required. 

In our opinion, a most interesting aspect of the IDM is the large flux of positrons, anti-protons, and anti- 
deuterons that are predicted to be produced in the annihilation of a light WIMP, with Mdm ^ 10 GeV. This 

^The simplest and most natural extension of the IDM is to add heavy singlet neutrinos with odd parity [8] but this model does 
not lead to charged lepton pair annihilation channels. A possibly interesting alternative might be to embed the IDM in a Universal 
Extra Dimension (UED) model (see for instance [6]), in which case annihilation in charged lepton pairs is allowed. 



14 



is due to two simple features. The first one is that the annihilation rate is large, dominated by Higgs exchange, 
in order to reach the WMAP abundance. The second feature is that the flux is proportional to the number 
density squared, which is large for a light candidate. One issue in confronting a light WIMP candidate to 
observations is that the low energy cosmic rays are very sensitive to solar modulation, a phenomenon which 
is still poorly understood. If progress could be made on this issue, our impression is that useful constraints 
might be put on the properties of light dark matter candidates, with Mbm ^ 10 GeV, a category that 
actually encompasses many models of dark matter. In the meantime, a prediction that is not obscured by 
solar modulation is that the flux of anti-deuterons produced by a light IDM dark matter candidate is predicted 
to be above the expected AMS-02 and GAPS experiments sensitivities . 
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